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Abstract 

We calculate the multi-instant on corrections to the ground state energy in large N Ma- 
trix Quantum Mechanics. We find that they can be obtained, through a non-perturbative 
difference equation, from the multi-instanton series in conventional Quantum Mechanics, 
as determined by the exact WKB method. We test our results by verifying that the one- 
instanton correction controls the large order behavior of the 1/N expansion in the quartic 
potential and in the c = 1 string. 
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1 Introduction 

The 1/N expansion of U(N) gauge theories plays a central role in our current under- 
standing of nonperturbative gauge theory dynamics. Most importantly, in some cases 
this expansion can be reinterpreted as a genus expansion in a dual string theory, leading 
to a connection between gauge theories and gravity theories. 

The 1/N expansion is in general an asymptotic expansion, and it has non-perturbative 
corrections of the form 0(e~ N ) due to large N multi-instantons. These corrections can 
trigger large N phase transitions [21, 12] and, by general arguments [16], they should 
control the large order behavior of the 1/N expansion. Moreover, when the gauge theory 
has a string dual, they can be reinterpreted in terms of D-branes. In spite of their 
relevance, the calculation of exponentially small corrections to the 1/N expansion has not 
been pursued in detail, even in exactly solvable, low dimensional models. 

The simplest toy models for the 1/N expansion are matrix models and Matrix Quan- 
tum Mechanics. Both models were solved at the planar level in the pioneering paper by 
Brezin, Itzkyson, Parisi and Zuber [5], and they have found an increasingly larger range of 
applications. Matrix models are now ubiquitous in physics, with important applications 
to two-dimensional gravity and topological string theory. However, explicit formulae for 
non-perturbative corrections in generic (off-criticality) matrix models have only appeared 
quite recently [14, 19, 18, 20, 22]. 

Matrix Quantum Mechanics has also been useful in many different areas. It provides 
for example a nonperturbative definition of c = 1 strings [10] and of one-dimensional 
superstrings [8, 23], and it has been instrumental in understanding some aspects of the 
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AdS/CFT correspondence [4]. The purpose of this paper is to present explicit formulae for 
multi-instanton corrections in Matrix Quantum Mechanics, by focusing on the energy of 
the ground state, which is the basic observable of the theory The first step in calculating 
these corrections is to find a useful characterization of the 1/N expansion to all orders. 
The next-to-leading correction to the planar result of [5] was found in [25], and it is not 
difficult to work out the full expansion. It turns out that, as in the case of matrix models 
[18], the ground state energy is determined by a difference equation. This equation relates 
the Matrix Quantum Mechanics problem, at all orders in the 1/N expansion, to the WKB 
expansion of the ground state energy in standard Quantum Mechanics, and for the same 
potential. We can then use the beautiful results obtained in nonperturbative Quantum 
Mechanics, specially in [26, 6, 27], in order to extract the full large N multi-instanton 
series, at all loops, in Matrix Quantum Mechanics. 

To illustrate our method, we analyze in detail Matrix Quantum Mechanics in a quartic 
potential, and we verify numerically that the one-instanton amplitude computed with our 
methods (at one loop) controls the large order behavior of the 1/N expansion. We also 
consider the double-scaling limit of this potential and reproduce in this way the large 
order behavior of the free energies of the c = 1 string. 

The paper is organized as follows. In section 2 we review the exact WKB method in 
quantum mechanics, including multi-instanton corrections, following mostly the results of 
[6]. In section 3 we start with a review of large N Matrix Quantum Mechanics, and we 
reformulate the 1/N expansion in terms of a difference equation. We then proceed to the 
determination of the multi-instanton corrections, and finally we analyze the connection to 
large order behavior. We conclude with some open problems. In the Appendix we include 
for completeness a short review of the nonperturbative treatment of the Schrodinger 
equation of [6]. 

2 Exact quantization conditions in Quantum Mechan- 
ics 

As we will see in the next section, the 1/N expansion in matrix quantum mechanics, 
as well as its non-perturbative corrections, can be calculated in a rather simple way by 
using the nonperturbative WKB method developed in [26, 24, 27, 6, 7]. In this section 
we review some relevant results from this method. 

2.1 Perturbative quantization conditions 

Let us consider the standard time-independent Schrodinger equation 

h 2 (p"(x) + p 2 (x, E)<p(x) = 0, p(x, E) = y/2(E - V[x)). (2.1) 
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If we write the wavefunction as 



(p(x) = exp 



Y(x')dx' 



(2.2) 



we transform the Schrodinger equation into a Riccati equation 

Y 2 (x)-ih^^-=p 2 (x,E), 
which we solve in power series in h: 

oo 

Y(x,E,h) = Y,Yk(x,E)h k . 

fc=0 

The functions Yk(x,E) can be computed recursively as 
Y (x,E)=p(x,E), 

Y n+1 (x,E) = 



If we split Y(x,E,h) into even and odd powers of h, 

Y(x, E, h) = F odd (x, E, h) + P(x, E, h 2 ), 

we find that 



Y odd (x,E,h) 



ih P'(x,E,h 2 ) 
~2 P(x, E, h 2 ) ' 



and the wavefunction reads 

ip(x, E, K) = 



y/P(x,E,V) 



i /- P(x',E,h 2 )dx' 



(2.3) 



(2.4) 



(2.5) 



(2.6) 
(2.7) 

(2.8) 



Let be a relative minimum of the potential V(x), and let us consider an energy E 
such that there are two roots of E — V (x) = near x*. These are classical turning points 
for the potential. Let B be a contour that encircles these roots, as depicted in Fig. 1. 
The classical frequency for the oscillation between these two turning points is 



ujb(E) = j> p(x,E)dx 
and the Bohr-Sommerfeld quantization condition reads 

uj b (E„) = 2nh(n - n = l,2, 
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(2.9) 



(2.10) 



V(x) 




Figure 1: A contour encircling the classical turning points. 

It will be convenient in the following to write 

X = h(n- 1 -) 

so that the solution to (2.10) is given by a function Eq°\x). 
Example 2.1. For the quartic potential 

V{\)= l 2 \ 2 + g\ i 
with g > 0, the Bohr-Sommerfeld frequency is given by 

/a 
d\ v / 2E - A 2 - 2#A 4 , 
■a 

where a, —a are the classical turning points given by 



^lQEg + 1 - 1 
49 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



The integral in (2.13) can be explicitly computed in terms of the elliptic functions E(k) 
K(k) as 

lu(E, g) = -(2g)^(a 2 + b 2 fAb 2 K(k) + (a 2 - b 2 )E(k) 



where 



, 2 y/16Eg + 1 + 1 



k 



Ag a 2 + b 2 ' 

The Bohr-Sommerfeld quantization condition leads to the energy levels 

E( X ) = X+ \gx 2 - ^9 2 X 3 + 0(g 4 ) 



(2.15) 
(2.16) 

(2.17) 



as a power series m g. 
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The condition (2.10) only incorporates the leading order WKB solution. An all-orders 
quantization condition was first proposed by Dunham in [9], and it reads as follows. We 
first define 

» oo 

n p (E,h)= i P(x,E,h 2 )dx = lo b (E) + V^ n) (E)n 2n , (2.18) 

Jb n=l 

where 

Jb\e) = £dxY 2n (x,E). 
The all-orders quantization condition reads 

and the solution is a power series in H, 



EM( X ,h)=Y, E{ nW n . 



(2.19) 



(2.20) 



(2.21) 



n=0 



The discrete energy levels are obtained by setting x to its quantized values (2.11). 

Example 2.2. In some cases, the quantization condition (2.20) makes possible to compute 
the exact energy levels [3]. Let us consider the potential 



V(x) = 



V 



cosh x 



(2.22) 



In this case, 



1 \n 



u B (E)=2ni K ^/2V -V^2Ej, uf{E) = 2n^/2V \^J (—J , (2.23) 
The infinite sum appearing in the l.h.s. of (2.20) can be summed up by using that 

oo 



n=0 



n) \8V 



2V + 



(2.24) 



and the quantization condition gives the exact energy levels 



E(x,h) 



X-(2^„4) 



2 X l 



(2.25) 



5 



2.2 Nonperturbative quantization conditions 



Nonperturbative corrections to energy levels in Quantum Mechanics can be obtained from 
a nonperturbative version of the quantization condition (2.20). This condition involves a 
function 

Q(E,h) = n p (E,h) + n np (E,h) (2.26) 

where the perturbative piece is given by (2.18), and Q np (E,h) is nonperturbative in h. 
The exact quantization condition reads 

Q(E( X ,h),h) = 2tt X . (2.27) 

and it results in nonperturbative corrections to the energy levels. Physically, fl np (E, H) is 
due to multi-instantons. 

The calculation of Q np (E,h) has been developed in a systematic and rigorous way 
in the context of the theory of resurgence [6], building on previous work (notably by J. 
Zinn- Justin) on multi-instantons in Quantum Mechanics [26, 27]. In the Appendix we 
summarize the approach of [6], which is quite elaborated. In this section we will illustrate 
the general method by considering standard situations where multi-instantons play a role: 
the case of an unstable potential (or false vacuum), and the case of a double- well potential. 

In general, Q(E,h) can be expressed in terms of the so-called Voros multipliers a 1 . 
A Voros multiplier is labelled by a contour 7 on the Riemann surface of the multivalued 
function of x, p(x, E). It is defined by 

P(x,£,/i 2 )dx}. (2.28) 
and we write 

00 „ 
P(x,E,h 2 ) = ^ + J2^ n) ( E ) h2n > ^\E) = f Y 2n (x,E)dx . (2.29) 

n=l T 

Let us first consider a false quantum-mechanical vacuum, like the one depicted in 
Fig. 2. In this kind of situation, as it is well-known, the energy levels develop an imaginary 
part which reflects the instability. This imaginary part is nonperturbative in H, and in 
order to compute it we have to compute the full Q(E, H) in (2.26). The perturbative part 
(the same as (2.18)) can be written in terms of Voros multipliers as 

tt p (E,h) = y loga B . (2.30) 

Let us assume that the potential has an instability associated to an extra turning point, 
as depicted in Fig. 2. Let us call A the cycle in the Riemann surface going from the 
turning point at the unstable minimum to the turning point of the instability. The 



a ' = exp 
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Figure 2: The A cycle goes from the turning point at the unstable minimum to the turning 
point of the instability. 



nonperturbative correction to fl(E, h) can be obtained by using the techniques developed 
in [6] and summarized in the Appendix. It can be written in terms of the Voros multiplier 
for the A-cycle as [6] 

Q np (E, h) = - { log(l + a A ) = -iha A (l - ^ + ^ + . . . V (2.31) 
At leading order, 

a A « e"^, u A = j ip{x)dx. (2.32) 
Notice that uja is real and positive for real E. 

Remark 2.3. As reviewed in the Appendix, the quantities Q(E, h) are defined by Borel 
resummations. In the case of an unstable potential, and it is well-known, the Borel 
transform has a singularity in the positive real axis and one has to use a prescription to 
avoid it. In the formulae above we use right Borel resummation as described in Appendix 
A. Thus Q in (2.27) is actually Q + , corresponding to the right Jost symbol J + (E) = 
1 + a B {l + a A ). The left Jost symbol is J~(E) = 1 + a B and is purely perturbative. 

The exact, nonperturbative quantization condition is then given by (2.27), where the 
perturbative and the nonperturbative part are given, respectively, by (2.30) and (2.31). 
Let us represent the solution of this condition as 

E( X ,h) = J2E (e \x,h), (2.33) 

where 



n -° (2 34) 

E^( X ,h) = he-^ E(0) y^ E n ] (xW- 



n>0 
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In this series, En\x) is the (n + l)-th loop correction in the £-th instanton sector. Eq°\x) 
is for example the result one obtains from the Bohr-Sommerfeld quantization condition. 
The full solution for the energy levels is then a double series in h, e _aM//?i , which is an 
example of trans-series solution. 

Let us calculate for example E^\x^f\) at leading order in h. After inserting (2.33) 
into the equation (2.27) and identifying the powers of exp(—u A /h), we obtain 

^|(4 0) (X))£ (1) (X, h) + n Qp (E^\ X )) + • • • = 0, (2.35) 
so that at leading order in K we find 

B o'(x,n) = a(^)" .-»"^«»(i + o(n)) 



T(Ef\ X )) 

where 



[I> e-W4 o) (x))(i + (ft)) 



(2.36) 



is the period of the classical trajectory of energy E between the turning points. 

Higher order corrections to (2.36) can be obtained recursively. In fact, it is possible 
to write down a compact and relatively explicit expression for the solution of (2.27). The 
function E(x) solves the equation 

uj b {E) + ij)(E) = 2n X (2.38) 

where 

oo 

${E) = uj b (E) + n Qp (E), cvl(E) = h 2n J B ] . (2.39) 

n=l 

Then one can write 

1 dE(x) 1 f dX f _±mjl d\ 



/cu b (A)+^(A)-2tt X /' 



2tt d X lo' b (E)+^'(E) J u B (X) +^(A) -2ir X J u B (\) - 2ir X 



_e_ wj?\x)) 1 -DME^ 

— ; e d x 2ir ; — ; e U V^ 



k=0 



(2.40) 

where on the second line E^ is the solution to the equation uj b (E) = 2nx (i.e. is the 
energy defined by the Bohr-Sommerfeld quantization condition), and 

D = ^ = nEf)dEf' (2 ' 41) 
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Here the normal ordering : : is denned by the rule "all £)'s to the left". One can now act 
by D^ 1 on both sides of (2.40) to obtain the explicit solution 

E(y) =:D~ 1 e -^ B o 0) ) : _J_ = : D -i e -^(4 0, )(i + a A ) mD : (2.42) 

One can extract from this solution the expression for the £-instanton correction: 

E d) = | P (£)( hD ) . e -°<4 : ^ (2.43) 



where 



Since cj^ = 0(h 2 ), and 



where 



'"w^ nM), l>l! (2 - 44) 

fc=i 



(fifl)e-^ = 5=e-^(l + O(fi)), (2.45) 
T,„„(B) = -^#> (2.46) 



BE 

we easily deduce the one- loop expression for a general multi-instanton correction £ > 1: 



° ~T(i#>)<! I r(25<°>) 1 ' ( } 










^ a 





Figure 3: The inverted quartic potential. The A cycle goes from a to b. 



Example 2.4. Let us consider the inverted quartic potential (2.12) with g 
this case, A = is a false vacuum, see Fig. 3. We first write 



-K < 0. In 



2E — X + 2/tA 4 = 2«(a a - \ 2 )(b 2 - A ) 



(2.48) 
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where 



2 -VI - IQEk + 1 2 VI - 16Ek+ 1 



a = — : , & = " : • (2-49) 

a, —a are the classical turning points, while b, —b are the turning points of the instability. 
The relevant period integrals are, 



u B = 2{2k) 1 / 2 j v^a 2 - A 2 )(6 2 - A 2 )dA, 

J —a 

oo A = 2{2k) 1 / 2 I V( A2 -a 2 W -A 2 )dA. 

J a 



(2.50) 



The periods can be computed with elliptic functions. The elliptic modulus is 

k 2 = (2-51) 



and we obtain 
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Ub = -(2«)36[(a 2 - b 2 )K(k) + (a 2 + b 2 )E(k)\ , 
uj a = \(2K)h\(a 2 + b 2 )E(k') - 2a 2 K(k')] , 



(2.52) 



where as usual k' 2 = 1 — k 2 . Of course, since the potential is symmetric, there is another 
contribution to Q np . It comes from the Voros multiplier associated to the cycle which 
goes from —a to —b. Therefore 

Sl(E, h) = j log a B + 2 ^ log(l + a A ) . (2.53) 

To check this result, obtained with the WKB approximation, we can compare it to the 
leading order result (in k) for the one-instanton correction to the energy of the n-th level. 
This is presented in e.g. [6, 13] and reads (we set H — 1) 

/ , \ n-l/2 

£«(«) = -— i (-) e-&(l + 0(K)) (2.54) 

n y ! V2^(n-1)! \kJ K V ,} y J 

To compare this with the result using the WKB method presented above, we have to keep 
all orders in h, but only the leading terms in k. We find that 



u>A = ±-X*>g% + 0( K ), iu#> = -£_ + 0(«), (2.55) 



where 

1 7 31 127 

Cl " 24 ' ° 2 ~ 2880 ' C3 " 40320 ' ° 4 " 215040 ' ' ' ' (2 ' 56) 
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The Voros multiplier a A is then given, in this limit, by 



and one finds 



aA (x, k) = exp | - i- - X log f e + ^ + 0{k) | , (2.57) 
E W (X^) = ^[1) e- Xl ° gf+ £^e-^(l + 0(K)). (2.58) 



One can check that 



-!°g(— Tf^J ="Xlog- + L^i, (2-59) 



as an asymptotic expansion in powers of x, an d thus expressions (2.58) and (2.54) coincide 
when we set x = n ~ 1/2- 




Figure 4: The double well. The A cycle goes from —a to a. 

It is possible to use the formalism of [6] to study other cases, like the double-well 
potential shown in Fig. 4. As it is well known (see for example [27]), the energy levels 
split into even and odd levels, according to the symmetry properties of the wavef unctions, 
which we will denote by E n ,±, respectively. Like before, the perturbative contribution to 
Q is given by (2.18), where the B cycle encircles one of the degenerate minima of the 
potential. Define now 

1 + ia A/2 

Q =wh (2 ' 60) 

where A is the cycle shown in Fig. 4. The exact energies E±(x,h) are defined by the 
quantization conditions 

n ± (E ± ( X ,h)) = 2n X , (2.61) 
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where 



Q±{E) = n p (E)±ih\ogQ. 



(2.62) 



This is a compact way of encoding the non-perturbative quantization conditions for the 
double- well potential first conjectured by Zinn-Justin [26]. One can check that the leading 
contribution to the energy difference coincides with the known answer 

E + -E- = ™e-'-P(l + 0(h)), (2.63) 

where 

s(E) = j \p(x,E)\dx (2.64) 

J —a 

and —a, a are the turning points for the A cycle. 

Remark 2.5. In this case, the results are obtained by using the so-called median resum- 
mation described in Appendix A. The result of this resummation is manifestly real. The 
computation of the left and right Jost symbols can be found in [6]. The median Jost 
symbol is 

j>°* = 2 + (l + a A )*(a B + a- B ) oc (1 + e^ n+ )(l + e^ n -) . (2.65) 



3 Nonperturbative effects in large N Matrix Quan- 
tum Mechanics 

3.1 1/N expansion of the ground state energy 

In Matrix Quantum Mechanics (MQM) the degrees of freedom are the entries of a Her- 
mitian N x N matrix M, and the Euclidean Lagrangian is 



J M 



Tr 



-M 2 - V N (M) 



where V N (M) is a potential. Notice that this problem has a symmetry 

M -> I MC 



(3.1) 



(3.2) 



where U is a constant unitary matrix. MQM can be regarded as a one-dimensional field 
theory for a quantum field M(t) taking values in the adjoint representation of U(N). 

As first shown in [5], the ground state energy of MQM has a 1/N expansion which 
can be obtained in terms of a system of free fermions. In this section we will review the 
results of [5] and we will extend them to all orders in the 1/N expansion. The Hamiltonian 
operator of MQM is given 



H = Tr 



1 d 2 

2dM 2 



+ V(M) 



(3.3) 
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where 

f) 2 f) 2 

Tr— = V — (3.4) 

OM 2 l^dM ab M ba 1 ; 

ab 

In order to study the spectrum of this Hamiltonian, it is useful to write the matrix M as 

M = UAU* (3.5) 

where 

A = diag(Ai, A 2 , • • • , \ N ) (3.6) 
is a diagonal matrix. It is easy to show that (see for example [1]) 

where 

A(A) = Y[(\ a - At) (3.8) 

a<b 

is the Vandermonde determinant, and T a b are differential operators w.r.t. the angular 
coordinates in U. 

Let us now consider singlet states. These are invariant under the U (N) group, and in 
particular they depend only on the eigenvalues A a up to permutation. The reason is that, 
after reduction to eigenvalues, the U (N) group still acts through the Weyl group, i.e. 
by permuting eigenvalues. Therefore, singlet states will be represented by a symmetric 
function, 

tt(Ai). (3.9) 

If we are now interested in computing the spectrum of the Hamiltonian for singlet states, 
we can reformulate the problem as a problem of iV fermions in the potential VJv(A). To 
see this, we introduce a completely antisymmetric wavefunction 



The equation 
can now be written as 

where h(X) is the Hamiltonian 



$(A) = A(A)*(A) (3.10) 

= (3.11) 

N . 

5>(A0 m\ j ) = E^(X j ) (3.12) 



h(X) = ~^ + V N (X). (3.13) 
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Since the fermions are not interacting, we can just solve the Schrondiger equation for a 
single particle of unit mass, 

/i(A)0 n (A) = E n <t> n {\). (3.14) 

In particular, the ground state of the system (in the singlet sector) will be obtained by 
putting the N fermions in the first N energy levels of the potential, and its energy will be 



N 



S(N) = J2En (3.15) 



71=1 



We want to compute the ground state energy S(N) at large N, and as an expansion 
in 1/JV. In order to have a good large N limit, Vn(\) must be of the form [25], 

V N (X) = NV(X/VN). (3.16) 

After rescaling 

A — > iV^A, (3.17) 

the Schrodinger problem becomes 

1 d 2 



2iV 2 dA 2 



+ V(X) U n (A) = e n n (A). (3.18) 



where we denoted 

e n = ^E n . (3.19) 

In this equation, 1/N plays the role of h, and this suggests using the WKB approximation 
in the calculation of the energy levels e n . The total energy of the ground state is now 

N N 
k=l k=l 

It will be convenient to re-scale the potential in such a way that the Schrodinger 
problem (3.18) becomes 

t 2 d 2 1 



2y2(U2 • r(\)jo„i\) I 2 o„iA). (3.21) 

The value of a depends on the potential, and t will be later identified with the 't Hooft 
parameter. 

Example 3.1. For the potential 

V(X) = y +tA 4 , (3.22) 
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we find the form (3.21) after setting A — > X/t, and a — 1, and 

v(X) = y + A 4 . 

For the potential (2.22) with V = 1/t 2 , we have a = 0. 
We now set 



(3.23) 



e fc = t 2 a e k 



and 



t = g s N. 

We will write the 1/JV expansion of the ground state energy as 

oo 
9=0 

The Schrodinger problem (3.18) becomes 



(3.24) 
(3.25) 

(3.26) 



9sj^ + v(X) \(j> n (X) = e n n (A). 



If we denote by 



X 



n — 
~~N~ 



z = g s {n - 1/2) = tx, 



(3.27) 



(3.28) 



we find that the all-orders perturbative WKB solution for the energy levels is given by 

n p (e(z,g 8 ),g a ) = 2Trz, (3.29) 
and it defines the perturbative function 



e(z,g s ) = ^2e n (z)g 2 s n . 



(3.30) 



?1=0 



For z = t, i.e. x = 1 or equivalently n ~ N, the function eo(t) is the WKB approximation 
to the Fermi energy of the fermionic system [5]. For future use, we will denote it by ep{t). 



Example 3.2. For the potential (2.22) we have 

1 



-V / 2+l4 



(3.31) 
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We now derive a general formula for £(t,g s ), which generalizes [5, 25] to all orders 
in the 1/N expansion. We first notice the following analogue of the Euler-Maclaurin 
asymptotic formula, 



N 



n=l 



n-1/2 
N 



N 



e n dx 0( x ) 



AT2 ( e d/dx _ !) _ 



1 

2N 



X =0 



sinh 



d/d X 
2N 



x=o 



N 2 / 0(x)dx-^AT 



2-2A; 



(1 - 2 1 " 2 *) fi 2fc d 2fe "V(x) 



k=i 



(2k)\ 



dx 



2fc-l 



X=l 



x=o 



Define first 

Then, one has that 

and using (3.32) one finds 



F(t,g s )=t 1 - a £(t,g s ). 



N 



N 

E 

n=l 



(3.32) 
(3.33) 

(3.34) 



^•^(^ = dz e(z, g s ) -J^g 



(l-2^ k )B 2k d 2k - 1 e(z,g s 



k=l 



(2k)\ dz 2 *- 1 
We can then write (3.35) as a difference equation 

9.{r{t + f)-r{t-f)} = e(t, g .). 



z=t 



z=0 



(3.35) 



(3.36) 



We then see that the large iV expansion of the ground state energy in MQM can be ob- 
tained, through (3.36), from the WKB expansion of the energies in an ordinary Quantum 
Mechanics problem with potential v(X). Moreover, (3.36) can be used as well to compute 
non-perturbative corrections, as we will see in a moment. 

In the case of a symmetric potential, like the double-well, the equation (3.36) has to 
be modified as follows: 



9s{r(t+^-r(t-?£)} = e + (t,g a ) + e4t,g a ), 
where e±(t,g s ) are obtained from the quantization condition (2.62). 



(3.37) 



Remark 3.3. The equation (3.36) is very similar, formally, to the equation determining 
the total free energy of the one-cut matrix model with the method of orthogonal polyno- 
mials, which can be also reformulated as a Toda-like difference equation [18]. The function 
e { z -,9s) plays the role of the function R(z,g s ), which is obtained as the continuum limit 
of the coefficients r n appearing in the recursion relation of orthogonal polynomials. 



16 



Example 3.4. For the potential V(A) in (2.22) we find, by a direct calculation using the 
result (2.25), 

£(t,g s ) = Jvf> = gfi-l - t + \\j<2^f] - ± (3.38) 

n=l * ' 

In this case, a = and 

F(t,g s )=t£{t,g s ) (3.39) 

satisfies indeed (3.36) with e(t,g s ) given by (3.31). For the quartic potential (3.22) we 
obtain from (2.17) the power series expansion for the planar approximation, 

which is the classic result of [5]. 

3.2 Nonperturbative corrections to the 1/N expansion 

The perturbative expansion of S(t, g s ) (or, equivalently, to J-(t, g s )) can be obtained from 
(3.36) by plugging in the perturbative expansion of e(t, g s ) (3.30). But the function e(t, g s ) 
has non-perturbative corrections coming from the exact WKB quantization condition 
(2.27). Therefore, we have 

e(t,g s ) = J2e (e) (t,g s ), (3.41) 

e>o 

with 

oo 

e^(t,g s ) = J2e 9 (t)g 2 s 9 , 



9=0 



(3.42) 



e W (t, 9s) = gjT***'™* eg>(t)fi, £ > 1, 

n=0 

just as in (2.34). In order to obtain the nonperturbative corrections to T(t, g s ), we simply 
have to consider the difference equation (3.36) as an exact statement, and plug in the full 
expansion (3.41). This leads to a multi-instanton expansion for J r (t,g s ) of the form 

T(t,g s ) = J2^ ie) ^9s) (3.43) 



£>0 



where 



9=0 



(3.44) 



jM{t,g a ) = e-^teVM'^rPWg?, £ > 1. 



n=0 
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To see how this works, let us calculate the one-loop, one-instanton correction ^^(t) 
in the unstable potential of Fig. 2. By plugging the expansion (3.44) in (3.36) we find, 
first of all, that the instanton action appearing in (3.44) is indeed UA(e F (t)). A simple 
calculation gives 

-2 S mh(i^).f?>( t )=e«"«, (3.45) 
and the right hand side can be read from (2.36). We then find, 

_ i de F (t) 



^ {t) = -J^-^T 77d — V' (3 - 46) 

^ dt sinh(j^) 

We can now use that 



dujA_dujA ^de F _ o du A /dE _ o J A dx/p(x,e F (t)) 
dt ~ dE [6F) dt ~ Z7T du B /dE ~ Zm f B dx/p(x,e F (t)) [ ] 

to finally write 

^ = 2mr{e F {t)\ (3.48) 

where t(E), given by the quotient of periods in (3.47), is the modulus of the curve defined 
by 

y 2 =p(x,E). (3.49) 
We then obtain an expression for ^^(t) of the form, 

^(t) = -(2||(e F (t))sin(vrr(e F (t)))) _1 , (3.50) 

which is written solely in terms of periods on the curve (3.49). In terms of the period of 
the trajectory T and the tunneling time 

T^- 9 ^ (3.51) 



we can also write 



2T M smh(^) 



(3.52) 



One can derive an expression for T analogous to (2.42). First of all, let us note that 
the relation (3.36) can be written as 

T = o • I ( 9 s a\ e (*> 9s) ■ (3.53) 

2<? s smh(f |) 

Then 

F = 9 n . \, m ■■ e- D <^-\l + a A (g s ,e F ))^ D : - 1 — , (3.54) 
2g s D smh{ng s D) T(e F ) 
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where 

1 d 

D = 7f^—^- (3.55) 
T(e F ) oe F 

The £-instanton contribution is 

{l) _ J_ P W(g a D) . , . {o-Y 
* -2&.SH&(*g a D) - 6 • T ' ^ 5bj 



with one-loop part 



^ - 2 T£!sinh(7r^) • (3 ' 57) 



3.3 Large order behavior of the I/TV expansion 

By standard arguments [16], the large order behavior of the 1/N expansion of the ground 
state energy should be governed by the instanton corrections that we have computed. 
Therefore, we can test our results by comparing them to the behavior of the amplitudes 
J- g {t), as g becomes large. We will restrict ourselves to the case of the unstable potential. 
The double-well potential, which is slightly subtler, can be studied in a similar way (see 
for example [27]). 

Let us write the one-instanton amplitude as 

oo 

g 2 s ^\t,g s ) = ig; b e-fsJ2^n(t)g:. (3.58) 

n=0 

Then, the perturbative amplitudes should have the large order behavior (see for example 
[19]) 

C—2g—b / ,, en— 1 \ 

In (3.58) the one-instanton amplitude is the discontinuity across the positive real axis, i.e. 
the difference between the results obtained with the right and the left Borel resummations. 
In our case, the left resummation is purely perturbative, and we can use the results for 
the one-instanton amplitude of section 3.2. We also have that b = —2, S = UA(ep(t)), 
and we obtain the leading asymptotics 

— y (2 r 2 > ( ^ M) :H a p.60) 

3 2tT(e F ) sinh(irT, un (e r )/T(e F )) 

We will now present test this formula with the inverted quartic potential (2.12) that 
we discussed above and its double-scaling limit. We will fix the normalization by choosing 
g = — 1. The energy is then related to the modulus k 2 introduced in (2.51) by 

£ « = (3 ' 61) 
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Er 



Figure 5: The critical point of Matrix Quantum Mechanics in the inverted quartic poten- 
tial takes place when the Fermi level 6f reaches the maximum E c . 



therefore the relation E(k) = ep(t) gives an implicit relationship between t, the 't Hooft 
parameter, and the modulus k. When k 2 = 1 we reach a critical point and the 1/N 
expansion breaks down. Physically, this critical point occurs when the Fermi level ep{t) 
attains the maximum of the potential E c , as shown in Fig. 5. This critical point plays 
an important role in non-critical string theory, since it makes possible to define the c = 1 
string by a double-scaling limit, see [10, 15] for reviews. 

We have computed ^^(t) up to g = 7 in order to test the large order behavior (3.60). 
Since the potential is symmetric, we have two identical contributions from instantons 
going from a to b, and from instantons going from —a to —b, therefore in this case we 
have to add an extra factor of 2 in (3.60). To do the test, we notice that the sequence 



UJE) =uT{E) sinh(7rT tun (E)/T(E)) 



(ME)) 29 ' 2 wo) 
(20-3)! ^ 



(3.62) 



should approach 1 as g — > oo, for all < k < 1. Of course we have only a few points to 
test this, but we can use Richardson transforms to accelerate the convergence. We recall 
that, given a sequence U g , g = 0, ■ ■ ■ , L, the r-th Richardson transform gives a sequence 
of L + 1 — r terms 

[R « u)n = E ^r*:': 1 '"" (3.63) 

m>0 



m!(r — m)\ 



with accelerated convergence. In Fig. 6 we plot the first few functions U g (k) as well as 
(R^U) 2 (k). The agreement is quite good, and confirms our analytic derivation of the 
one-instanton correction. 
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k 

U 2 (k) U 3 (k) [/„{*) U 5 (k) U 6 (k) — U 7 (k) 1 ( R (S) v )^ k) 



Figure 6: In this figure we plot the first functions U g (k) and the second element of the 
fifth Richardson transform (R^U)2(k). 

3.4 Application to the c = 1 string 

As we mentioned above, if we consider the inverted double- well potential with k — >• 1, we 
can define a double-scaled theory by considering the limit 

E c — e F = [i — > , — = z fixed, (3.64) 

This limit defines the c = 1 string, and the function !F^(t, g s ) becomes (see, for example, 
[15]) 

Therefore, in this limit the asymptotics as g — > oc can be computed directly by using that 

~ f^' 9 » 1- (3-66) 
The genus g coefficient of (3.65), which we will denote by J 7 ^, behaves like 



9 

1 

2~7 3 



J? ~ — (2n)- 2 ^T(2g - 2) (3.67) 
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Let us compare this to the predictions of (3.60). The different quantities involved in that 
expression can be easily computed in the double-scaling limit. Since 1 

E C -V(X) = -^(X-K) 2 + --- (3.68) 
Denoting £ = A — A* one computes 

T tun (e F ) = 4 j ^ = =2tt, 



o 



(3.69) 



-uj A (e F ) = - [ d^y^J^e = ^~ = — ■ 

7a 9s J 9s z 



9. 



Since T(ep) diverges in this limit, 

T(e F ) sinh^vr^M^) -> 7rT tun (e F ). (3.70) 

After including the extra factor of 2 to account for the symmetry of the potential (which 
is also included in the formula of [15]), we find that (3.60) reproduces (3.67). We should 
mention that already in the paper [11], the leading contribution to the large order behavior 
of the c = 1 string is determined by computing the action of the instanton in the double- 
scaling limit. The paper [2] also discusses nonperturbative effects in the c = 1 string. 



4 Conclusions and open problems 

In this paper we have determined the non-perturbative corrections to the ground state 
energy in large N Matrix Quantum Mechanics. Essentially, we reduced the problem to 
the calculation of multi-instanton corrections in conventional Quantum Mechanics, which 
can be obtained in turn through exact quantization conditions. 

There are various possible generalizations of our work. We have restricted ourselves 
to the ground state energy, but the nonperturbative corrections to energies of excited 
singlet states or non-singlet states (like the adjoint state analyzed in [17]) could be also 
analyzed with our techniques. It would be also very interesting to analyze the large 
N multi-instantons once fermionic degrees of freedom have been introduced, as in [1]. 
Finally, we have worked out in detail the example of the unstable quartic potential and 
its double scaling limit, and we have given all the necessary ingredients to understand the 
double-well potential. It would be interesting to work out this example in more detail, 
since is relevant to the understanding of type OB strings [23, 8]. 

1 Here we actually use a different normalization of the potential (2.12) with g = — |. 
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A Resurgence and nonperturbative quantization con- 
ditions 

In this Appendix we briefly review the nonperturbative treatment of the Schrodinger 
equation based on the theory of resurgence and described in detail in [6, 7]. 

In this approach, the WKB expansions satisfying the Schrodinger equation, such as 
(2.8), as well as other relevant functions written in terms of series in h, are regarded as 
so-called resurgent symbols, i.e. formal sums of the form 

v = J2^ e ~ f ( A - x ) 

where (p u are formal series in h: 



n=0 



The formal power series <p w (h) must satisfy the condition that their Borel transforms 

in— 1 

(n~l)\ 



71=1 



have only finite number of singularities in the positive real direction, so that the inverse 
transform (which is called a resummation of the series (fu{h)) can be defined by 

a + [ e-£&,(Od£ (A.4) 



where 7 is some contour starting from zero and going along the positive real direction 
which avoids these singularities 2 . There are two natural choices of 7: the contour 7_ 
avoiding the singularities from the left, and the contour 7+ avoiding the singularities from 
the right. The corresponding results give us the left and right resummations: S-ip^^H) 
and s + (p u (h). One can construct from them the left and the right resummations of the 
resurgent symbol ip: 



^2s-<p u e « , s + <p = ^2s+<p u e » (A.5) 



2 



If nothing is known about the growth of <pu(£) at infinity, the function (A.4) can be defined only up 



to a function /(ft) of hyperexponential decrease, i.e. 3C : |/| < Ce r 
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These resummation operators s + and s_ actually define isomorphisms of the algebra of 
resurgent symbols into the algebra of so-called extended resurgent functions. The linear 
space of WKB symbols will be denoted by WKB. The two-dimensional linear space of 
solutions to the Schrodinger equation, regarded as extended resurgent functions, will be 
denoted by S. Notice that the same resurgent function can be obtained from two different 
symbols by using the different resummations s±. One can define the action of the Stokes 
automorphism S on WKB by requiring the commutativity of the following diagram: 

WKB >WKB . (A.6) 

S 

It turns out that WKB-symbols are well defined only in the so-called Stokes regions 
of the complex x-plane. We will denote the space of WKB symbols in the Stokes region 
R by WKB(fi). The complex plane is divided into Stokes regions by Stokes lines. Stokes 
lines are lines starting or ending (or both) at critical points (zeroes of the momentum 
p(x,E)) along which e^ p(^ e )<^ decreases/increases fastest. Inside a Stokes region the 
space of WKB-symbols can be decomposed into the direct sum with respect to the choice 
of sign of the momentum (±p(x, E)): 

WKB(_R) = WKB (+p) (i?) © WKB ( - p) {R). (A.7) 

The Stokes automorphism respects this decomposition inside the Stokes region R. 

Consider now two different Stokes regions R and R' and choose some resummation 
prescription (s + or s_). Then there is a map 

C% RI : WKB(i?) -> WKB(#) (A.8) 

called the connection isomorphism, which relates two different WKB-symbols correspond- 
ing to the same global function through the different resummation s± . 

Let us now consider a simple pattern of Stokes lines: three unbounded Stokes lines 
L, L', L" going out from one critical point: 




(A.9) 



One can decompose 

WKB(i2) = WKB (+p) (i?) © WKB ( - p) (i?) (A. 10) 




24 



and 

WKB(-R') = WKB (+p) (i?') © WKB^^i?), (A. 11) 

where WKB* +!) ' consists of symbols dominant along L, i.e. with the factor 

exp(i J X p(£,E)d^ (A.12) 

increasing along L, and WKB^ P ' consists of symbols recessive along L, i.e. with the 
factor 

exp(~ J'p(£,E!)dz) (A.13) 

decreasing along L. Then the elementary connection isomorphism can be written in the 
following matrix form: 

C% R , = ( I I ) , (A.14) 

where L is analytic continuation across L and 5l is analytic continuation along the contour 
around L, as shown in (A. 9). The elementary connections operators for the right and left 
resummations are the same. 

The global picture of Stokes lines can not be always represented as a composition of 
simple patterns as the one depicted in (A. 9). In general one has bounded Stokes lines. 
However, one can always make small deformations such that all Stokes lines will be un- 
bounded, and one can then find connection isomorphisms C RR , between every two Stokes 
regions R, R' by using the elementary connection isomorphisms that we just described. 
Deformed Stokes lines are lines along which 

exp^ J X p(t,E)d^ (A.15) 

decreases/increases fastest, where e — > 0, see Fig. 7. One can use the deformation with 
e > (respectively, e < 0) to compute the connection isomorphism for the right (resp., 
left) resummation. 

The problem of finding energy levels can be formulated as finding the values of the 
energy E such that the subspace S-^ C S of solutions descending at +oo (or having 
negative momentum if one searches for resonances in unstable potentials) coincides with 
the subspace of solutions S +tXl C S descending at +oo (having negative momentum in 
the case of resonances). This is equivalent to the vanishing of the so-called Jost operator 
J : 5_oo — > S/S +00 . If one chooses the basis (Jost basis) (0, (jf) of S such that <fi E S-oo 
and 0* ^ S +OQ , the action of the Jost operator is given by a Jost function J(E) defined 
by 

4> = J{E)4>* mod 5 +0O . (A. 16) 

Then the energy levels are given by the equation J{E) = 0. 
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e<0 



e>0 



Figure 7: Stokes lines for the cubic potential and their resolution. 

Let us consider the WKB-symbol (p defined in the Stokes region Ri unbounded in the 
negative real direction, such that <fi is the resummation of tp. If 0* is the resummation 
of the WKB-symbol (p defined in the Stokes region R2 unbounded in the positive real 
direction, we can write 

Cr 1 r^ = J ± (E)0 mod s^S+oo. (A.17) 

J^E) is called the right/left resurgent symbol of the Jost function. Correspondingly, 
the right/left symbols of the energy levels are the roots of the equation J ± {E) = 0. It is 
usually convenient to represent Jost symbols in the form 

J ± = l + e* n± V S) . (A.18) 

The nonperturbative quantization condition is then 

n ± (E)=27v X . (A. 19) 

Thus the problem of finding a nonperturbative equation for the energy reduces to finding 
the connection isomorphism C^ iFt2 using the rule (A. 9). In many examples, the physically 
relevant resummation prescription is the median resummation defined by 

8^ = 8.00-5=8+003. (A.20) 

The corresponding median Jost symbol is 

jmcd = J- = g-i (A.21) 

This prescription has the property that, in the case of a stable potential, the solutions to 
the equation 

J™ d (E) = 1 + eR nmod = (A.22) 
are resurgent symbols with real coefficients. 
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